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Abstract 

We use a restricted sample of elliptical galaxies, whose kinematical 
parameters inside the semimajor axis were calculated correcting the 
effect of the integration of the light along the line of sight, in order 
to analyze a possible relationship between the mass of a Supermassive 
Black Hole (SMBH) and the kinetic energy of random motions in the 
host galaxy. We find M B h oc (M G a 2 ) a with 0.87 < a < 1 depending 
on the different fitting methods and samples used. This result could 
be interpreted as a new fundamental relationship or as a new way to 
explain the old Mbh — c law. In fact, the relations of the velocity 
dispersion both with the mass of the SMBH (Mbh oc er 4 ' 12 ) and with 
the mass of the host galaxy (Mq oc a 216 ) induce us to infer an almost 
direct proportionality: Mbh oc Mqo 2 . A similar relationship is found 
for the total kinetic energy involving the rotation velocity too. 



1 Introduction 



The discovery in the nucleus of an increasing number of galaxies of SMBHs 
and the estimation of their masses with several methods have allowed the 
beginning of statistical analyses that try to show the influence of a SMBH 
on the properties of the host galaxy. In particular, several relationships were 
proposed between the mass of a central SMBH and the velocity dispersion 
H] • [2] > EH > the bulge luminosity or mass [1] [3] or the dark matter halo [5] , of 
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the corresponding galaxy. Among them, the relationship with the smallest 
scatter is Mbh ot <J a where 3.75 < a < 5.3. The different values of a, found 
by several authors using different samples and different fitting methods, are 
well discussed in the paper of Tremaine et al. [3| that obtain 

LogM BH = (4.02 ± 0.32)Lo 5 (o-/200) + (8.13 ± 0.06) (1) 

where (and also afterwards in the rest of our paper) Mbh is expressed in 
solar masses, a in Km/s and logarithms are base 10. 

While the astronomers believe in these relationships, in such a way that 
they use them for example to predict the value of the mass of a SMBH in a 
galaxy of a distant cluster [7j, from the theoretical point of view there are 
many possible interpretations about the physical meaning and the origin of 
these relations jS] jH]. In this paper we verify the existence of a relationship 
between the rest energy of the SMBH and the kinetic energy of random 
motions inside the semimajor axis of the host elliptical galaxy Mbhc 2 oc 
{M G a 2 f. This could be a new fundamental relation or it can be used (in 
a simple way if (3 — 1) to interpret and explain the Mbh oc a a law (where, 
for example, a ~ 4 as in eq. (1), but, in principle, we could also find a 
different value) , through a hidden dependence of the mass of the galaxy on 

——2 

the velocity dispersion. In that case a relation Mq oc a? holds. Finally we 
consider the role of the rotation velocity in a relationship between the rest 
energy of the SMBH and the total kinetic energy. The paper is organized as 
follows: in sect. 2 we describe the sample of galaxies chosen for the statistical 
analysis and the model used to obtain the kinematical data, in sect. 3 the 
results of the fitting procedure are given and in sect. 4 the conclusions are 
drawn. 

2 The sample 

As suggested by Gebhardt et al. :2| , a particular care has to be devoted, in 
the choice of the sample, to find suitable kinematical parameters that must 
be referred not only to the central regions but for example to the effective 
radius, in order to have a higher signal to noise ratio and to be less sensitive 
to the direct influence of the gravitational field of the central SMBH. In 
our case, we want to consider not only the velocity dispersions but also 
the masses of the host galaxies. There are several methods to compute the 
masses of galaxies, for example using in different ways the Virial Theorem 
or starting from the Jeans equation or better with self consistent models 
and recently using also data on globular clusters, planetary nebulae and the 
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temperature and density distribution of the thermal X-ray gas. In the past, 
one of us (A.F.) used two of these methods. When the mass of the whole 
galaxy was to be considered, the Virial Theorem was used in its scalar form 
2T + U = and, knowing the rotation velocity and velocity dispersion, the 
corresponding Virial Mass was derived ^Uj • On the contrary, if only the part 
of the galaxy inside a radius R e is studied, the system is not really isolated 
and the effect of the matter at a radius r > R e has to be considered as 
an unknown external pressure jllj . This term and the rotation velocity are 
often neglected in the case of elliptical galaxies. In order to overcome this 
problem, Busarello et al.^2j preferred to compute the masses of a sample of 
62 elliptical galaxies inside the effective semimajor axis A e (related to the 
effective radius R e by the formula 

T 2 = 1 " ~ ( 2 ) 

where e is the ellipticity) , starting from the Jeans equation describing the 
equilibrium of a spheroidally symmetric system having eccentricity e and an 
isotropic velocity dispersion tensor: 

1 d . W y- (r)] - J^W _ - 4jG "- e2)1/2 r ^Mr, , (3) 



p(r) dr r r J (r 2 — x 2 e 2 ) 1 / 2 

where r is the radius in the equatorial {z = 0) plane, p(r) is the (unknown) 
spatial density, a(r) and v(r) are the one-dimensional velocity dispersion 
and rotation velocity respectively, and G is the constant of gravity |13| . 

The solution of Eq. (3) can be written in the form p(r) = poXl(r), where 
l(r) is the luminosity density. Busarello et al. |14j (hereafter BLF) assumed 
that the luminosity distribution corresponds to the spatial deprojection of 
the r 1 / 4 law. A simple analytical approximation for the deprojection of the 
r 1//4 law has been derived by Mellier and Mathez |15j : 

l(r) = r-Pexpi-br 1 ^) (4) 

where (3 = 0.855 and b = 7.669. Substituting this solution together with 
V(r) and a(r) in the equation (3) an expression of po as a function of r 
can be obtained. Computing the value of po for each object in the sample 
at 10 different radii, the residuals with respect to its mean value < po > 
turn out to be very small (and will be used to estimate the error AM), and 
show no systematic trend with the radius, thus supporting the hypothesis 
that po = constant at least in the inner regions (r < A e ) 12 . So, the final 
result for the mass density is p(r) =< po > r~^exp{—br l / A ). Sometime it 
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happens that two ways to estimate the same quantity give different results; 
actually, in our case we observe that this method to compute masses can 
lead to values that differ even 30% from masses obtained using the Virial 
Theorem. 

We decide to adopt along this paper all the kinematical parameters com- 
puted by Busarello et al. and published for a and V of 54 elliptical galaxies 
in BLF and for the mass and the specific angular momentum in ^H] • These 
results have the advantage to be treated with the same method and the 
same fitting procedure; furthermore they are all referred to the effective 
semimajor axis, as we required, and are corrected for one projection effect. 
The observable quantities are actually affected by two types of projection 
effects: the inclination of the rotation axis with respect to the line of sight 
and the integration of the light along the line of sight. Busarello et al. |14j 
|17j applied in fact a method to correct the second effect, deprojecting the 
rotation velocity curves and the velocity dispersion profiles. They start from 
a very simple model assuming that an elliptical galaxy 

1. has a spheroidal symmetry, 

2. follows the de Vaucouleurs r 1 / 4 law whose spatial deprojection was 
given by the approximated analytical expression (4), 

3. the rotation axis is perpendicular to the line of sight, 

4. the rotation velocity has cylindrical symmetry and 

5. the velocity dispersion is isotropic and has spherical symmetry. 

The method has been explained in details in BLF where also the sources 
of kinematical data (except for N3379 ^Sj and IC1459 ^1) an d the values 
of the parameters used in the calculations are listed. It allows to compute, 
starting from a fit of the experimental data, the deprojected V(r) and a(r) 
to be inserted in the equation (3) in order to obtain the masses. From 
those two analytic functions it is easy to obtain also the specific kinetic 
energies due to rotation and random motions respectively and some others 
kinematical parameters such as specific angular momentum, spin etc. 

In this paper we consider as the reference sample, the intersection be- 
tween the set of elliptical galaxies studied by Busarello et al. HH 
and the set of SMBH masses in Table 1 in the paper of Tremaine et al. 
The masses and the kinematical parameters of the resulting 14 galaxies are 
listed in table 1 where V is the luminosity weighted mean rotation velocity 
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inside A e (related to the rotational kinetic energy by Ty = MV 2 /2) and 



a 2 ^< a 2 yy >= 



St 



j y (r)pd 3 x = 
/ pd 3 x 3M 



(5) 



is the luminosity weighted mean of the line of sight component, of the veloc- 
ity dispersion tensor, assuming that the mass-to-light ratio is constant inside 
A e and that the tensor is isotropic (T a is the corresponding kinetic energy). 
We do not insert in our statistics the galaxy N221 because it has a mass 
two orders of magnitude less than all the others and probably a different 
dynamical behavior and its velocity dispersion is less than 70Km/sec so it 
belongs to the set of galaxies that even with modern observations can be 
neglected for problems of instrumental resolution |20| . Finally our sample 
is restricted to only 13 elliptical galaxies and, even if it seems small and 
derived from old data (but the sample of is smaller and the kinemati- 
cal data are as old as ours) and from a procedure "extremely sensitive to 
errors or incompleteness in the data" [2^ > its values of a give results in full 
agreement with the relationship (1). 



As a first test we plot the data in the [Log{a / '200), LogMsn) plane and 
the line in Fig. 1 is obtained by using a standard least-squares fitting and 
assuming, just like 0, "that errors in dispersion measurements are zero and 
that errors in LogMsn are the same for each galaxy" 1 . The best-fit line is 



that confirms the eq. (1) obtained by Tremaine et al.jHj- 

The linear correlation coefficient is r = 0.827. Furthermore we can 
calculate the unknown constant error in LogMsH using the formula: 



for a relation of the form y = a+f3x, and we obtain e y = 0.46. If we compute 
the average of the experimental errors in LogM^H (table 1 in [3] ) , we obtain 
< e y >= 0.19, so either there is an underestimation of the experimental 
errors (due for example to our hypothesis that neglects errors on a), or the 
intrinsic dispersion of the relation is larger than the measurement errors. 

1 These hypotheses will be relaxed later. 



3 The results 



LogM BH = (4.12 ± 0.84)Loc/(o-/200) + (8.35 ± 0.13) 



(6) 




i=l 
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TABLE 1. - The sample with the data used in the paper. 





(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Galaxy 


Type 


Mbh 


a; 


a 


V 


Ma 


AM S 






(m.) 




(km ■ sec ' 


(km ■ sec 1 


t«„> 


~McT 


N221 


E2 


2,5 ■ 10" 


43,9 


60 


37 


4,0 ■ 10 s 


0,07 


N821 


E4 


3,7 ■ 10' 


52,8 


180 


117 


1,4 • 10" 


0,04 


N2778 


E2 


1,4 ■ 10' 


18,3 


140 


96 


4,0 • 10"' 


0,16 


N3379 


E1 


1,0 ■ 10 s 


37,1 


193 


53 


3,0 ■ 10'° 


0,12 


N4291 


E2 


3,1 ■ 10 s 


15,9 


250 


76 


9,0 • 10"' 


0,05 


N4473 


E5 


1,1 ■ 10* 


31,4 


191 


62 


7,0 • 10"' 


0,16 


N4486 


EO 


3,0 ■ 10" 


108,9 


269 


20 


3,5 ■ 10" 


0,10 


N4564 


E3 


5,6 ■ 10' 


29,8 


125 


143 


5,0 • 10"' 


0,16 


N4649 


E1 


2,0 ■ 10' 


83,7 


224 


56 


2,2 ■ 10" 


0,51 


N4697 


E4 


1,7 ■ 10" 


92,5 


177 


151 


2,9 ■ 10" 


0,19 


N4742 


E4 


1,4 ■ 10' 


12,8 


81 


91 


1,0 ■ 10"' 


0,17 


N5845 


E3 


2,4 ■ 10 s 


4,6 


236 


73 


2,0 • 10"' 


0,16 


N6251 


E2 


5,3 ■ 10* 


10,0 


288 


54 


2,3 • 10" 


0,12 


IC1459 


E3 


2,5 ■ 10" 


46,3 


282 


383 


3,5 ■ 10" 


0,03 



NOTE. - The sources of the data are: columns (1) and (2) from [3], columns (3), (4) and (5) 
from BLF and columns (6) and (7) from [12, 16]. 
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Now our aim is to study the dependence of the mass of the galaxy inside 
the effective semimajor axis, on the corresponding velocity dispersion. From 
the dependence of the mass-to-light ratios on the luminosity M/L cx L ' 25 
and from the Faber - Jackson relation L oc a 4 , Ferrarese and Merritt 
infer M cx <t 5 for early type galaxies. On the other side, Busarello et al. [21] 
starting from the data derived in BLF (the same used in this paper) for a 
sample of 40 galaxies showed that "no real correlation holds between M/L 
and V when the masses are computed through the Jeans equation. Then, 
following the reasoning of Burkert and Silk we could expect from the 
application of the Virial theorem (Mq oc R e a 2 ) and from the results of the 
Fundamental Plane of elliptical galaxies (R e cx a 2 ) that Mq cx <t 4 , but the 
best-fit of our data gives: 

LogM G = (2.16 ± 0.71) Log (a/ '200) + (10.98 ± 0.11) (8) 

This relationship, even if it has a large scatter and a poor correlation (r = 
0.678), induces to think that 

M BH oc M G a 2 cx cj 4 - 16 (9) 

in agreement with the previous fit (6). 

In order to check this hypothesis, we show in Fig. 2 the best-fit of the 
relation 

Mr a 2 

Log(M BH ) = (0.87 ± 0.17)Log— %— + (4.33 ± 0.80) (10) 

& 

with a correlation coefficient r = 0.833 even better (of course with our data) 
than the famous relation (6). 

We have also tried to consider the role of the total kinetic energy adding 
the contribution due to the rotation velocity and we obtain: 

Mr (3a 2 + V 2 ) 

Log(M BH ) = (0.88 ± 0.18)Log^± — ^— + (3.80 ± 0.92) (11) 

c z 

with a correlation coefficient r = 0.829 and the fit is interesting too (Fig. 3). 
It is remarkable that the slope of the relation (11) is the same as the relation 
(10) and that the intercept changes in the way we expect. Furthermore, we 
can calculate the unknown constant error in LogM B n using (7) and we ob- 
tain again e y = 0.46 for both the relations (10) and (11) confirming that the 
procedure leads to a value 2.4 times greater than the average experimental 
error. So we must check what happens if we use a fit weighted by the ex- 
perimental errors and take into account also the errors on the independent 
variable. 
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As for our kinematical data it is not available the estimate of the error on 
a for each galaxy, we have considered until now the optimistic assumption of 
a negligible error. Now it is time to check if something changes in the relation 
(10), considering the worst scenario of a relative error 20% on a as we can 
estimate from the discussion on the possible sources of error (anisotropy, 
triaxiality, fitting procedure, inclination, etc.) contained in BLF. To this 
aim we use the effective variance method suggested by Orear [221 that with 
an iterative procedure minimize the x 2 given by the formula 



X 



2 _ (Vi-a- [3xi) 



^(Ay,) 2 + /? 2 (Ax,) 2 



for a relation of the form y = a + (3x. The corresponding results are listed in 
Appendix A and show no change in the slope of the relationships found until 
now, the only advantage being a decrease of the difference between the resid- 
uals and the corresponding experimental errors. On the contrary, slightly 
different slopes are obtained starting from the Akritas and Bershadi [22] 
statistical method (used also in pP) with the corresponding 68% confidence 
interval: 

Mr a 2 

Log(M BH ) = (0.98 ± 0m)Log— %— + (3.8 ± 0.4) (12) 

tr 

with a reduced x 2 ( see Appendix A) x 2 = 1>9. 

This fit is shown in fig. 2 with a dashed line and involves Mbhc 2 oc 
Mqo~ 2 and could be considered as interesting as the famous relation (1) is. 

The same check for the relation (11) gives 

Mr (3a 2 + V 2 ) 

Log(M BH ) = (0.995 ± 0.082)Log— ^ — ^— '- + (3.21 ± 0.42) (13) 

cr 

(dashed line in fig. 3) with x 2 = 1-8. 

Even better results can be obtained using a reduced sample (hereafter 
SAMPLE B) of galaxies if we eliminate the two ellipticals with the largest 
residuals N821 and N4697. The relationships obtained applying the fitting 
procedure to the remaining 11 galaxies are listed in Appendix B. The sat- 
isfying results are the increase of all correlation coefficients, the decrease of 
all the values of Xr an d the slopes of the two relationships that involve the 
kinetic energies that are closer to the unity. While a coefficient less than one 
is more difficult to be interpreted, it would be more interesting from the the- 
oretical point of view, the meaningful possibility of a direct proportionality 
between the mass of a SMBH and the kinetic energy of random motions of 
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the host galaxy. Moreover, if our results are confirmed, the relation between 
Mbh and Mq will be non linear [H] |25| . 

The standard scenario of a Black Hole life predicts that a part of the total 
mass of a SMBH is due to the accretion process and a part of this last mass 
can be converted in radiation and ejected again in the region surrounding 
the Hole. Prom our results it seems that a part of the rest energy of the 
Black Hole is in this way strictly related to the kinetic energy of random 
motions of the stars of the host galaxy even in a region far from the Hole. 

4 Conclusions 

From our sample of data we derive the suggestion to consider a relationship 
between the masses of SMBHs and the Kinetic Energy of random motions, 
or even the total kinetic energy, in elliptical galaxies. Only a deeper analysis 
with the sophisticated machinery today available and with new data |26| . 
could discriminate if our interpretation can work and if it is generally valid or 
it holds only for a restricted sample of galaxies (for example the ellipticals 
with the mass in a small range of values). Given the assumptions of the 
model, it is clear that triaxiality, anisotropy of the velocity field, inclination, 
deviations from the r 1//4 law, are all sources of possible errors that could 
affect the derived kinematical parameters. In BLF some of these errors 
were discussed in detail and their total influence on the final results was 
estimated in no more than 20%. However the underestimation of one of 
the above effects, the restricted sample and even the way to compute the 
masses we used, can lead to draw wrong conclusions. On the other side we 
think that it would be surely worse to neglect a priori the suggestion that 
comes to us by the relationships from (10) to (13) and above all from (10B) 
to (13B). So we have now arguments to ask again the question contained in 
the title: is there a relationship between the mass of SMBH and the kinetic 
energy in its host elliptical galaxy? 
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Appendix A 

The formula used in this paper to estimate all the maximal errors in the 
functions F of the parameters (a,b,c, ....) shown in table 1, is 



AF(a,b, c) 



dF 



da 



Aa + 



dF 



db 



Ab + 



dF 



dc 



Ac 



but our results do not change if we use for example the formula of Ferrarese 
and Merritt (cited also by Tremaine et al. 3 ) to estimate the error bar 
in LoqMbh that is 

ower) 

IXLoglviBH = t. 
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In order to compare our results with the ones obtained by [I] , we have used 
in (12) and (13) the same fitting method |23| and the reduced x 2 given by 
the formula: 

2 _ 1 \ - (yi-a- (3xj) 2 
Xr N - I ^ (Ay^ + P^(A Xi y 

for a relation of the form y = a + (3x. On the other side, the results obtained 
using the iterative procedure |22] (the errors Ax, are ignored in the first step) 
are the following: 

LogM BH = (4.10 ± 0.73)Log(a/200) + (8.41 ± 0.12) (6A) 

with Xr = 11; 

LogM G = (2.15 ± 0.38) Log {a/ 200) + (10.97 ± 0.06) (8A) 
with Xr = 3.0; 

Mr a 2 

Log(M BH ) = (0.88 ± 0.06)Log— %— + (4.30 ± 0.46) (10A) 
with Xr = 2.2; 

Mr (3a 2 + V 2 ) 

Log(M BH ) = (0.89 ± 0.07)Log^± — ^— + (3.77 ± 0.52) (11A) 

c z 

with Xr = 2.1. 



Appendix B 

We neglect the two ellipticals with the largest residuals: N821 and N4697. 
The relationships contained in the main part of the paper are reanalyzed 
using the remaining 11 galaxies forming the so called SAMPLE B. The 
results obtained using the iterative procedure [22] are the following: 



LogM BH = (4.06 ± 0.73)Log(o-/200) + (8.46 ± 0.13) (65) 
with a linear correlation coefficient r = 0.846 and Xr = 0.9; 

LogM G = (2.26 ± 0.52) Log (a/ 200) + (10.90 ± 0.09) (8B) 
with r = 0.755 and Xr = 1-6; 

Mro~ 2 

Log(M BH ) = (0.92 ± 0.07)Lo#— ^ + (4.25 ± 0.49) (10B) 
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with r = 0.907 and xl = 0.9; 

Mr (3a 2 + V 2 ) 

Log(M BH ) = (0.92 ± 0.07)Log— ^ — ^- '- + (3.70 ± 0.55) (115) 

cr 

with a r = 0.907 and xl = 0.9. 

The results obtained using the Akritas and Bershadi method are the 
following: 

Mr a 2 

Log(M BH ) = (1.03 ± Qm)Log—^ + (3.71 ± 0.38) (125) 

c l 

with xl = 1-06; 

s , s M g (3ct 2 + T/ 2 ) , , 
Log(M BH ) = (1-06 ± 0.08 )Log—^ — = + (3.02 ± 0.39 135 

c z 

with xl = 1-01. 
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FIG. 1 



Figure 1: SMBH mass versus luminosity weighted, line of sight velocity 
dispersion of the host elliptical galaxy. The solid line is the best linear fit 
with the least squares method (6). The error bars for Mbh are from [3] 
and for a are calculated from the upper limit on the relative error of 20% 
estimated from BLF. 



FIG. 2 



Figure 2: SMBH mass versus Mqo 1 /c 2 of the host elliptical galaxy. The 
solid line is the best linear fit (10) with the least squares method, the dashed 
line is the best fit (12) with the Akritas - Bershadi method [23]. The error 
bars for Mbh are from for Mq are in table 1 and for a are calculated 
from the upper limit on the relative error of 20% estimated from BLF. 



FIG. 3 



Figure 3: SMBH mass versus Mg(3cf 2 + V 2 )/c 2 of the host elliptical galaxy. 
The solid line is the best linear fit (11) with the least squares method, the 
dashed line is the best fit (13) with the Akritas - Bershadi method |23| . The 
error bars for Mbh are from for Mq are in table 1 and for a and V are 
calculated from the upper limit on the relative error of 20% estimated from 
BLF. 



